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[i]  The  first  flyby  of  Mercury  by  the  Mercury  Surface,  Space  Environment, 

Geochemistry  and  Ranging  (MESSENGER)  spacecraft  occurred  on  14  January  2008. 

In  order  to  provide  contextual  information  about  the  solar  wind  (SW)  properties  and  the 
interplanetary  magnetic  field  near  the  planet,  we  have  used  an  empirical  modeling 
technique  combined  with  a  numerical  physics-based  SW  model.  The  Wang-Sheeley-Arge 
(WSA)  method  uses  solar  photospheric  magnetic  field  observations  (from  Earth-based 
instruments)  in  order  to  estimate  inner  heliospheric  conditions  out  to  21.5  solar  radii 
from  the  Sun.  This  information  is  then  used  as  input  to  the  global  numerical 
magnetohydrodynamic  model,  ENLIL,  which  calculates  SW  velocity,  density, 
temperature,  and  magnetic  field  strength  and  polarity  throughout  the  inner  heliosphere. 

The  present  paper  shows  WSA-ENLIL  conditions  computed  for  the  several  week  period 
encompassing  the  first  flyby.  This  information  is  used  in  conjunction  with  MESSENGER 
magnetometer  data  (and  the  only  limited  available  MESSENGER  SW  plasma  data)  to 
help  understand  the  Mercury  flyby  results.  The  in  situ  spacecraft  data,  in  turn,  can  also  be 
used  iteratively  to  improve  the  model  accuracy  for  inner  heliospheric  “space  weather" 
purposes.  Looking  to  the  future,  we  discuss  how  with  such  modeling  we  can  estimate 
relatively  continuously  the  SW  properties  near  Mercury  and  at  the  cruise  location  of 
MESSENGER  now,  for  upcoming  flybys,  and  toward  the  time  of  spacecraft  orbit 
insertion  in  2011. 
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1.  Introduction 

[2]  The  Mercury  Surface,  Space  Environment,  Geochem¬ 
istry,  and  Ranging  (MESSENGER)  spacecraft  was  launched 
on  3  August  2004.  The  spacecraft  then  began  its  long, 
circuitous  route  toward  eventual  insertion  into  Mercury 
orbit  in  March  2011.  MESSENGER  had  two  flybys  of 
Venus  (October  2006  and  June  2007)  and,  on  14  January 
2008,  had  its  first  of  three  flybys  of  Mercury.  This  successful 
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initial  transit  through  the  near  environs  of  Mercury  [Solomon 
et  al. ,  2008]  gave  the  first  close  views  of  the  Mercury  surface, 
atmosphere,  and  magnetosphere  since  the  Mariner  10  flybys 
of  1974  and  1975  [e.g.,  Russell  et  al,  1988]. 

[3]  During  the  period  around  MESSENGER'S  closest 
approach  to  the  planet,  which  occurred  at  ~1905  UTC  on 
14  January,  the  magnetometer  and  the  Energetic  Particle  and 
Plasma  Spectrometer  (EPPS)  recorded  a  variety  of  interest¬ 
ing  magnetospherie  phenomena  [Slavin  et  al,  2008].  How¬ 
ever,  the  in  situ  measurements,  as  well  as  the  inferences 
from  measurements  made  prior  to  and  following  the  mag¬ 
netospherie  passage,  suggested  that  the  solar  wind  (SW)  and 
interplanetary  magnetic  field  (IMF)  B  during  the  particular 
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Figure  1.  Elements  of  the  real-time  coupled  eoronal- 
heliospherie  model  used  in  the  present  study.  GONG  data 
arc  used  as  input  for  the  synoptic  maps,  and  forecast  outputs 
are  given  at  the  test  bed  site  (http://Helios.swpe.noaa.gov/ 
enlil/latest-veloeity.html).  Definitions  and  abbreviations  are 
described  in  the  text. 

time  of  the  flyby  were  not  very  conducive  to  strong 
magnetospherie  substorm  or  particle  acceleration  events  of 
the  type  seen  by  Mariner  10  instruments  [e.g..  Baker  et  ai^ 
1 986]  due  to  the  fact  the  IMF  was  northward  throughout  the 
passage. 

[4]  The  limitations  of  having  measurements  of  the  inner 
heliospheric  SW and  IMF  conditions  at  only  a  single  location 
have  made  it  evident  that  a  broader  contextual  view  of 
the  Mercury  (and  MESSENGER)  space  environment  is 
valuable.  For  example,  it  is  useful  to  know  about  large- 
scale  and  medium-scale  interplanetary  structures  that  may 
be  present  near  MESSENGER,  such  as  high-speed  SW 
streams,  eorotating  interaction  regions,  and  IMF  polarity 
reversal  boundaries.  Present-day  solar  observations  (made 
from  Earth-based  facilities),  along  with  empirical  and 
physies-based  models,  have  the  potential  to  provide  a  good 
global  picture  of  the  inner  heliosphere  [Arge  et  al.^  2004; 
Owens  et  al.,  2005]. 

[5]  In  this  paper,  we  first  describe  current  state-of-the-art 
modeling  techniques.  We  then  show  the  application  of  these 
techniques  to  the  MESSENGER  flyby  period  in  order  to 
place  the  spacecraft  results  into  a  broad  physical  perspec¬ 
tive.  We  make  specific  comparisons  of  the  model  results 
with  measurements  at  1  AU  and  also  from  the  relevant 
MESSENGER  sensors  before  and  after  the  planetary  en¬ 
counter.  We  conclude  with  a  general  discussion  of  a  strategy 
for  using  our  modeling  approaches  for  future  MESSENGER 
and  other  (remote)  Mercury  observations. 

2.  Modeling  Inner  Heliospheric  Properties 

[6]  The  Wang-Sheeley-Arge  (WSA)  model  is  a  combined 
empirical  and  physies-based  representation  of  the  quasi¬ 
steady  global  SW  flow  [Arge  et  ai,  2004].  It  has  been  used 
in  prior  work  to  predict  the  ambient  SW  speed  and  IMF 
polarity  at  Earth  (as  well  as  other  points  in  the  inner 
heliosphere)  and  is  an  extension  of  the  original  Wang  and 
Shceley  model  [Wcifig  and  Sheeley,  1992].  The  model  uses 


ground-based  observations  of  the  Sun’s  surface  magnetic 
field  as  input  to  a  magnetostatic  potential  field  source 
surface  model  [Schatten  et  al,  1969]  of  the  coronal  field 
(see  Figure  1).  The  effects  of  outward  flows  in  the  corona, 
which  are  not  explicitly  contained  in  the  formulation,  are 
approximated  by  the  imposition  of  radial  field  boundary 
conditions  at  the  source  surface,  which  is  a  Sun-eentered 
sphere  typically  positioned  at  2.5  Rs^  where  Rs  is  the  solar 
radius.  A  number  of  important  changes  have  been  made  to 
the  Wang-Sheeley  model  [Arge  and  Pizzo,  2000],  including 
(1)  improvements  to  the  model  inputs,  which  are  the  line  of 
sight  (LOS)  photospherie  magnetic  fields,  and  (2)  improve¬ 
ments  to  the  empirical  kinematic  model  itself 

[7]  The  photospherie  field  observations  are  the  basic 
large-scale  observable  used  to  drive  the  computations  and 
serve  as  a  key  input  to  all  coronal  and  SW  models.  Here  we 
use  updated  photospherie  field  synoptic  maps  (i.e.,  synoptic 
maps  updated  four  times  a  day  with  new  magnetograms) 
constructed  with  LOS  magnetograms  from  the  National 
Solar  Observatory’s  Global  Oscillation  Network  Group 
(GONG)  system.  These  ground-based  data  provide  the 
present  “standard”  data  input  for  WSA  modeling.  Because 
observational  evidence  suggests  that  the  solar  magnetic  field 
is  nearly  radial  at  the  photosphere  (except  in  strong  active 
regions),  the  LOS  field  measurements  from  these  data 
sources  have  been  converted  to  radial  [see  Arge  et  al.^ 
2004]. 

[r]  Using  WSA  results  relatively  near  the  Sun,  an  ideal 
magnetohydrodynamie  simulation  called  ENLIL  [Odstrcil 
et  al,  2004;  Toth  and  Odstrcil,  1996]  is  then  performed  to 
model  the  SW  flow  outward  to  beyond  1  AU.  The  compu¬ 
tational  domain  is  the  uniform  grid  occupying  the  sector  of  a 
sphere  defined  by  pairs  of  boundaries  at  fixed  radii  (inner 
and  outer),  at  fixed  meridional  angles  (north  and  south),  and 
at  fixed  azimuthal  angles  (east  and  west).  The  position  of 
the  inner  boundary  is  set  at  0.1  AU  (^^^21.5  Rs).  and  the 
outer  boundary  is  set  at  l.l  AU.  The  meridional  and 
azimuthal  extents  span  30  150°  and  0  360°,  respectively. 
The  inner  boundary  lies  in  the  supersonic  flow  region,  near 
the  outer  field  of  view  of  the  Large  Angle  and  Speetrometrie 
Coronagraph  C3  on  the  Solar  Heliospheric  Observatory 
spacecraft.  The  outer  boundary  is  chosen  to  enable  record¬ 
ing  of  simulated  temporal  profiles  of  SW  properties  at  and 
near  the  Earth  position  [Odstrcil  et  ai,  2004]. 

3.  Mercury  Flyby:  Modeled  Inner  Heliosphere 
Conditions 

[9]  The  product  of  the  combined  WSA-ENLIL  modeling 
is  a  specification  of  the  SW  flow  speed,  plasma  density,  SW 
plasma  temperature,  and  magnetic  field  strength  throughout 
the  inner  part  of  the  heliosphere.  Prior  modeling  work  has 
tended  to  focus  on  optimization  of  modeled  results  near  the 
Earth’s  location  at  1  AU,  or  at  the  first  Lagrangian  point, 
LI,  which  is  the  location  of  the  Advanced  Composition 
Explorer  (ACE)  spacecraft.  New  data  now  available  from 
the  dual  spacecraft  Solar  Terrestrial  Relations  Observatory 
(STEREO)  mission  have  provided  a  very  useful  broader  range 
of  SW  measurements  with  which  to  compare  WSA-ENLIL 
results  (see  http://stereo.gsfe.nasa.gov). 

[10]  A  color  representation  of  the  radial  flow  speed  (F^)  in 
the  equatorial  plane  computed  for  the  entire  inner  heliosphere 
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Figure  2.  (a)  Modeled  radial  SW  speed  and  (b)  density  in  the  equatorial  plane,  viewed  from  the  north 
ecliptic  pole,  obtained  from  the  WSA-ENLIL  model  near  the  time  of  the  first  MESSENGER  flyby  of 
Mercury.  The  color  scale  for  is  given  by  the  color  bar  above  Figure  2a.  As  shown  at  the  bottom,  the 
locations  of  Earth,  STEREO-A,  STEREO-B,  Venus,  Mercury,  and  the  MESSENGER  spacecraft  arc  all 
indicated  by  small  colored  dots.  The  inner  domain  of  the  model  (where  WSA  is  utilized)  is  denoted  by  the 
white  central  circle.  The  computational  domain  of  the  ENLIL  simulation  is  shown  by  the  colored  area.  A 
nominal  computed  magnetic  field  line  is  indicated  by  the  heavy  dashed  line  extending  from  the  inner 
boundary  outward  through  the  Earth  past  1  AU.  The  red-blue  color  coding  along  the  edge  of  the  outer 
boundary  of  computation  shows  the  polarity  of  the  IMF:  red  indicates  IMF  positive,  or  pointing  away 
from  the  Sun,  while  blue  indicates  negative  polarity  with  the  IMF  pointing  toward  the  Sun.  The  white 
curves  mark  the  estimated  IMF  polarity  sector  boundaries  in  the  equatorial  plane.  The  SW  density  in  the 
inner  heliosphere  modeled  by  WSA-ENLIL  and  shown  in  Figure  2b  is  scaled  by  to  the  value  at  1  AU. 


early  on  14  January  2008  is  shown  in  Figure  2a.  From  the 
model  results,  we  see  that  a  well-developed  broad  SW 
stream  region  was  expected  in  the  heliosphere  during  this 
time.  The  prominent  modeled  SW  speed  enhancement  (up 
to  ~600  km/s)  at  ~1  AU  was  in  the  approximate  longitude 
seetor  -135  0°  (mostly  behind  the  azimuthal  location  of 
the  Earth),  where  O'^  is  the  direction  toward  Earth.  According 
to  the  model,  the  stream  would  have  enveloped  STEREO-B 
at  the  time  of  the  snapshot  but  had  not  yet  reached  STEREO- 
A  (see  section  3.1  and  Figures  3  5  for  detailed  model- 
data  comparison).  MESSENGER  and  Mercury,  of  course, 
were  effectively  collocated  at  the  time  of  the  snapshot  and 
were  subjected  to  essentially  identical  SW  flow  conditions. 
From  a  Mercury  magnetospheric  response  perspective,  there 
was  no  significant  SW  speed  enhancement  expected  to  be 
seen  near  the  planet  on  the  particular  day  of  the  spacecraft 
flyby.  The  high-speed  stream  under  discussion  would  have 
passed  over  the  planet  (and  MESSENGER)  several  days  prior 
to  the  flyby  of  Mercury  with  the  highest-speed  (600  km/s) 
stream  features  having  been  expected  to  rotate  over  the 
Mercury  location  some  10  11  days  earlier.  Benign  SW 


conditions  at  Mercury  on  ~14  January  (as  suggested  by  the 
model)  meant  that  the  magnetosphere  was  relatively  quies¬ 
cent  during  the  spacecraft  passage,  as  was  observed  [Slavin  et 
aL,  2008]. 

[11]  Another  important  control  on  Mercury’s  magneto- 
spheric  configuration  and  physical  extent  is  the  density  (N) 
of  the  incident  SW.  Higher  density  along  with  higher  flow 
speeds  would  imply  higher  SW  dynamic  pressure  (F'dyn  “ 
mNV^,  where  m  is  the  effective  SW  ion  mass  and  V  is  the 
scalar  flow  speed)  and  also  could  suggest  higher  impacts  of 
SW  particles  onto  the  magnetic  polar  regions  of  Mercury’s 
surface  [sec  McClintock  et  ai,  2008].  The  WSA-ENLIL 
model  is  able  to  provide  global  modeling  estimates  of  SW 
density  throughout  the  inner  heliosphere  for  the  same 
modeling  domain  (as  was  shown  above  for  SW  speed). 

[12]  Color-coded  density  computation  results  for  14 
January  2008  are  shown  in  Figure  2b.  From  Figure  2b,  it 
is  seen  that  two  relatively  high-density  regions  of  SW 
plasma  were  expected  from  this  model.  One  lay  ahead 
(between  0  and  --45°  ecliptic  longitude  at  1  AU)  of  the 
high-speed  SW  stream  shown  in  Figure  2a  and  bracketed 
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Figure  3,  ACE  measurements  of  (first  panel)  SW  speed,  (second  panel)  density,  and  (third  panel) 
plasma  temperature  and  (fourth  panel)  IMF  magnitude,  all  as  1-h  averages  measured  at  the  upstream 
Lagrangian  (LI)  point,  arc  shown  by  the  red  traces  for  1  31  January  2008,  The  smooth  blue  curves  show 
WSA-ENLIL  model  values  for  the  same  period.  The  48  h  bracketing  the  first  MESSENGER  flyby  of 
Mercury  (vertical  line)  arc  shaded. 


LMJL-2.5  mcdfcs  VV^A-1.C  GONG  2008^01  “01  OOiOOlOO  -'•outun  ul  pu'^omrlm  ot  SILHLO^A 


J]-Cy^  01-05  C'l-OB  ai-U  01-17  01-21  iM-25  01 -2a  G2-U2 

TIME 


Figure  4.  Similar  to  Figure  3  but  for  the  STEREO-A  spacecraft  near  1  AU  at  solar  longitude  ahead  of 
the  Earth’s  position  (as  described  in  the  text). 
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Figure  5.  Similar  to  Figure  3  but  for  the  STEREO-B  spaeeerafl  near  I  AU  at  solar  longitude  behind  the 
Earth’s  position  (as  deseribed  in  the  text). 


the  IMF  polarity  reversal  boundary  (shown  by  the  white 
spiral  eurve  in  Figure  2b).  This  density  enhaneement  was 
modeled  to  be  greatest  near  1  AU,  but  it  extended  all  the  way 
inward  to  the  inner  boundary  of  the  ENLIL  eomputation. 

[13]  A  seeond  modeled  high-density  region  is  seen  in 
Figure  2b  to  braeket  the  other  IMF  polarity  reversal  region 
in  the  longitude  seetor  from  135°  to  ~150°  at  1  AU.  The 
boundary  of  a  third  slightly  inereased  density  region  was 
seen  in  the  model  ealeulations  to  extend  over  the  loeation  of 
Mereury  and  MESSENGER  at  the  time  of  this  snapshot. 
However,  this  density  stream  was  mueh  weaker  than  the 
other  two  that  were  in  other  seetors  of  the  inner  heliosphere. 

4.  Comparison  of  Model  Results  With  Data 
at  1  AU 

[14]  As  noted  in  seetion  3,  WSA-ENLIL  modeling  results 
have  been  used  extensively  to  “foreeast”  the  SW  and  IMF 
values  at  I  AU  using  ground-based  solar  input  data  (e.g., 
from  the  GONG  network).  In  prineiple,  sueh  modeling 
results  ean  give  3  4  day  foreeasts  of  SW  properties  [e.g., 
Baker  et  al.,  2004].  The  model  results  ean  be  readily 
eompared  with  real-time  measurements  from  sueh  spaee- 
eraft  missions  as  ACE  and  STEREO. 

[15]  The  period  of  January  2008  was  eharaeterized  by  two 
broad  SW  streams  (see  Figure  2)  that  were  persistent  and 
well  developed  throughout  the  inner  heliosphere  (Figure  3). 
ACE  observations  for  that  period  show  the  onset  of  a  high¬ 
speed  stream  on  5  January,  with  a  large-density  spike  just  at 
the  leading  edge  of  the  fast  stream.  Another  high-speed 
stream  eommeneed  on  about  13  January  with  a  weaker 
density  enhaneement  at  the  leading  edge  of  this  stream  as 
well.  This  seeond  stream  persisted  broadly  for  over  a  week 


at  ACE,  and  the  SW  speed  eventually  diminished  to  about 
400  km/s  by  ~23  January. 

[16]  The  WSA-ENLIL  modeling  results  for  this  interval 
of  time  are  shown  by  the  smooth  blue  lines  in  eaeh  panel  of 
Figure  3.  It  is  seen  that  most  of  the  general  features  of  the 
SW  properties  were  seen  similarly  in  the  model.  Certainly 
the  two  SW  streams  were  elearly  simulated  by  the  model. 
The  seeond  stream  (13  21  January)  is  very  well  fit  both  in 
time  profile  and  in  aetual  SW  speed.  During  this  period,  the 
density,  temperature,  and  magnetie  field  strength  values  all 
are  also  matehed  well  by  the  model.  The  first  stream  (5 
II  January)  is  less  well  fit  than  the  seeond  stream.  It  is 
seen  that  the  modeled  SW  speed  rises  too  late  (by  about 
1  day),  and  the  other  parameters  (N,  T,  and  B)  are  also 
modeled  to  rise  a  bit  late,  and  absolute  values  are  not  fit  as 
aeeurately  as  they  are  for  the  second  stream  event. 

[17]  It  is  interesting  to  compare  and  eontrast  model  results 
with  data  from  other  platforms  near  I  AU  for  this  same 
period.  We  show  in  Figures  4  and  5  the  measurements  at 
STEREO-A  and  STEREO-B  that  are  comparable  to  the 
measurements  from  ACE  in  Figure  3.  (The  relative  posi¬ 
tions  of  the  STEREO  spaeeeraft  are  shown  in  Figure  2). 
Because  of  the  general  physical  proximity  of  the  STEREO-A 
and  STEREO-B  spaeeeraft  to  ACE,  we  sec  most  of  the  SW 
stream  features  at  all  three  spacecraft.  However,  substantial 
timing  differences  arc  seen  at  the  three  locations. 

[is]  We  emphasize  that  correct  characterization  of  the 
solar  wind  stream  properties  at  the  three  separated  spatial 
locations  (ACE,  STEREO-A,  and  STEREO-B)  is  an  impor¬ 
tant  validation  of  the  modcFs  capabilities.  Getting  the 
arrival  times  correct  at  the  three  disparate  solar  longitudes 
indicates  that  the  solar  wind  stream,  indeed,  has  the  essen¬ 
tial  structure  predicted  by  the  model  and  shown  in  Figure  2. 
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Figure  6.  Computed  SW  parameters  for  January  2008.  The  shading  in  eaeh  panel  shows  the  interval 
immediately  surrounding  the  Mercury  encounter  period  for  MESSENGER  (13  15  January).  The  first, 
second,  and  third  panels  show  the  SW  speed,  density,  and  temperature  for  this  period,  respectively  (which 
matches  the  time  period  of  Figures  3  and  4).  The  fourth  panel  shows  the  ealeulated  value  of  the  IMF 
magnitude  for  the  same  interval  of  time.  The  red  data  points  plotted  in  the  first,  third,  and  fourth  panels 
show  MESSENGER  magnetie  field  and  available  plasma  data  (as  deseribed  in  the  text)  for  the  period 
1  January  to  1  February  2008  for  eomparison.  The  bottom  three  panels  show  other  (derived)  quantities. 
The  fifth  panel  shows  the  sonic  and  Alfven  Maeh  numbers.  The  sixth  panel  shows  estimates  of  the  eross- 
magnetospherie  potential  £drop  (see  text),  and  the  bottom  panel  shows  dynamic  pressure  P^yn- 


As  was  the  ease  for  the  ACE  comparisons  with  the  model 
output  (Figure  3),  WSA-ENLIL  does  a  good  job  of  fore¬ 
casting  the  SW  and  IMF  data  both  ahead  and  behind  the 
Earth  locations  as  seen  by  the  STEREO-A  and  -B  spacecraft 
comparisons. 

[19]  We  note  in  Figures  3  5  that  there  are  challenges  in 
getting  a  proper  eomparison  with  the  measured  temper¬ 
atures.  For  example,  on  STEREO  the  temperature  plots 
are  from  an  instrument  that  measures  only  protons  (not 
electrons).  Protons  have  higher  temperatures  than  electrons 
(especially  in  fast  streams  from  coronal  holes).  Numerical 
simulations  provide  the  “mean”  plasma  temperature  (which 


by  the  above  arguments,  should  therefore  be  lower  than  the 
plotted  proton  temperatures). 

5.  Comparison  of  Model  Results  With 
MESSENGER  Data 

[20]  The  SW  temporal  profiles  and  IMF  parameter  values 
that  were  ealeulated  from  the  WSA-ENLIL  model  at  the 
MESSENGER  location  from  1  January  to  1  February  2008 
are  shown  in  Figure  6.  From  the  model  results,  we  see  that 
the  time  of  MESSENGER  closest  approach  to  the  planet 
was,  indeed,  expected  to  be  a  period  of  relatively  slow  SW 
flow  and  relatively  high  density.  The  magnetie  field  strength 


6  of  8 


AlOlOl 


BAKER  ET  AL.:  SPACE  ENVIRONMENT  OF  MERCURY 


AlOlOl 


was  computed  to  be  gradually  varying  throughout  this  entire 
period,  but  the  IMF  was  at  a  local  plateau  (5  15  nT) 

during  the  encounter  period. 

[21]  In  order  to  compare  model  results  with  the  in  situ 
spacecraft  measurements,  we  plot  the  nearly  continuous 
MESSENGER  magnetometer  data  [Anderson  et  al,  2008] 
on  the  model  curve  of  B  in  Figure  6.  These  comparisons 
show  that  the  model  values  are  in  good  agreement  with  the 
measured  values  from  MESSENGER  throughout  most  of 
the  month-long  interval.  An  exception  was  for  a  brief  period 
on  20  22  January,  when  the  MESSENGER  field  magnitude 
dropped  dramatically  and  model  results  were  very  different 
from  the  measurements.  However,  in  broad  terms,  the 
ENLIL-MESSENGER  magnetic  field  agreement  was  quite 
good. 

[22]  Because  of  the  mounting  position  of  the  EPPS 
plasma  analyzer  [Zurbuchen  et  ai,  2008]  on  the  spacecraft 
and  because  of  spacecraft  pointing  constraints,  we  do  not 
have  continuous  or  complete  MESSENGER  measurements 
of  the  solar  wind.  During  limited  intervals  throughout  the 
period  covered  in  Figure  6,  however,  the  team  was  able  to 
obtain  a  sufficient  fraction  of  the  distribution  function  to 
generate  estimates  of  SW  parameters  (particularly  speed  and 
temperature).  Data  for  such  intervals  arc  plotted  in  red  on 
the  respective  model  curves  in  Figure  6.  The  MESSENGER 
speed  estimates  are  higher  than  the  ENLIL  calculations,  but 
the  inferred  temperatures  are  low  compared  to  the  model. 
Naturally,  the  measured  SW  parameters  show  more  struc¬ 
ture  and  higher  time  variability  than  the  model,  which 
would  be  expected  given  the  slower  cadence  of  inputs 
(ground  based)  and  inherent  spatial  smoothing  to  the  model 
calculations.  (Interestingly,  most  complete  plasma  measure¬ 
ments  were  derivable  during  the  time  of  lowest  and  most 
disparate  magnetic  field  values). 

[23]  The  fifth,  sixth,  and  seventh  panels  of  Figure  6 
provide  derived  parameters  of  considerable  relevance 
(and  utility)  for  magnetospheric  modeling  of  the  Mercury 
system.  The  fifth  panel  shows  the  sonic  Mach  number  (in 
blue)  and  the  Alfven  Mach  number  (in  green).  These  values 
(<10  throughout  the  encounter  interval)  are  pertinent  to 
estimating  the  expected  bow  shock  and  magnetopause 
properties  at  Mercury  [Slavin  et  al.,  2008].  The  sixth  panel 
shows  our  calculated  values  of  the  potential  drop  ^drop  (in 
kilovolts)  across  the  Mercury  magnetosphere.  We  have  used 
the  ENLIL  model  SW  speed  combined  with  the  measured 
magnetic  field  normal  component  and  a  scaling  distance  of 
1  Mercury  radius  (R^  =  2349  km)  to  compute  E^rop-  Such 
potential  drop  estimates  arc  an  important  aspect  of  the  kinds 
of  particle  acceleration  that  might  occur  within  the  Mercury 
magnetotail  [e.g.,  Baker  et  al,  1986].  The  seventh  panel 
shows  the  SW  dynamic  pressure  Pdyn- 

6.  Discussion  and  Future  Work 

[24]  The  combined  WSA  and  ENLIL  models  provide 
contextual  perspective  on  the  SW  and  IMF  conditions  in 
the  inner  heliosphere  during  the  time  around  the  first 
MESSENGER  flyby  of  Mercury.  Our  comparisons  both 
before  and  after  the  time  of  closest  approach  show  that  the 
MESSENGER  magnetometer  measurements  were  in  re¬ 
markably  good  agreement  with  model  results  for  most  of 
the  4-weck-long  interval  of  comparison.  One  relatively  brief 


period  (20  22  January  2008)  showed  completely  different 
magnetic  field  strengths  than  expected  from  the  model 
predictions.  Because  of  poor  sensor  placement,  available 
MESSENGER  SW  measurements  are  of  very  limited  utility 
for  comparison  with  model  outputs.  The  model  results 
clearly  show  that  Mercury’s  magnetosphere  on  14  January 
2008  was  being  subjected  to  an  extensive  region  of  rela¬ 
tively  low-speed,  quiet  SW.  This  information  helps  explain 
the  benign,  inactive  magnetosphere  that  was  sampled  by 
MESSENGER  sensors  [Slavin  et  ai,  2008]. 

[25]  The  magnetic  field  data  during  the  20  22  January 
period  showed  a  large  reduction  (generally)  in  the  radial 
(from  the  Sun)  component  and  also  more  fluctuations  in  the 
transverse  components.  As  seen  in  Figure  6,  the  overall 
magnitude  of  the  field  was  about  half  of  the  value  expected 
from  the  model  results.  Since  the  MESSENGER  plasma 
analyzer  was  also  able  to  measure  more  of  the  SW  plasma 
distribution  during  much  of  this  interval  (see  Figure  6),  wc 
infer  that  the  direction  of  SW  flow  must  have  been  different 
(and  more  propitious  for  measuring  its  properties)  during 
much  of  the  anomalous  magnetic  field  interval. 

[26]  Wc  have  looked  at  solar  images  and  other  remote 
sensing  data  for  the  period  around  20  January  2008  to  sec  if 
any  features  on  the  Sun  might  explain  the  weak  field 
interval  seen  in  the  MESSENGER  data.  It  is  not  obvious 
that  any  coronal  or  solar  wind  stream  feature  can  be 
discerned.  Thus,  this  aspect  of  the  data-model  comparison 
remains  a  puzzle. 

[27]  We  assert  that  the  present  results  are  important  in 
several  respects.  First,  the  methods  utilized  here  can  be  used 
in  a  similar  (or  improved)  fashion  for  future  MESSENGER 
flybys  of  Mercury.  Thus,  wc  will  be  able  to  supply  daily 
updated  values  of  forecasted  context  information  (~  1  2  days 
in  advance  at  Mercury)  for  future  in  situ  measurements.  This 
will  help  in  the  prompt  analysis  and  interpretation  of  flyby 
data  and  for  the  orbital  phase  of  the  MESSENGER  mission. 
A  second  point  is  that  the  MESSENGER  IMF  (and  any 
available  SW)  data  in  the  inner  heliosphere  arc  quite  useful 
local  “ground  truth”  for  the  WSA-ENLIL  model  calcula¬ 
tions.  Having  actual  observations  at  ^^.0.3  0.4  AU  helio¬ 
centric  distance  that  can  be  “assimilated”  into  the  model 
will  help  to  improve  and  guide  the  overall  model  perfor¬ 
mance.  in  this  sense,  data-theory  closure  can  lead  to  a  better 
overall  space  weather  prediction  at  Earth  by  the  WSA- 
ENLIL  combination  [e.g..  Baker  et  al,  2004]. 

[28]  As  a  final  point,  wc  can  look  forward  to  future 
measurements  made  of  the  Mercury  system.  For  example. 
Earth-based  measurements  of  the  atmosphere  of  Mercury 
[e.g.,  Killen  et  al,  2004]  often  would  benefit  from  having  a 
general  knowledge  of  what  type  or  intensity  of  SW  is  striking 
the  planet  (and  its  magnetosphere)  at  the  time  that  ground 
telescopic  data  arc  acquired.  The  modeling  shown  here  can 
help  provide  that  useful  contextual  information.  Beginning  in 
2011  when  MESSENGER  is  in  orbit  around  Mercury,  the 
spacecraft  will  be  within  the  magnetosphere  and  magnetotail 
of  the  planet  for  extended  portions  of  each  orbit.  Model 
results  will  provide  continuous  information  about  the  SW 
and  IMF  that  is  influencing  magnetospheric  dynamics  and 
exospheric  variability.  We  can  also  envision  using  the  WSA- 
ENLIL  time-dependent  specifications  and  forecasts  of  SW 
parameters  and  IMF  as  inputs  to  many  future  magnetospheric 
simulation  models  [e.g.,  Kahin  et  al,  2000]. 
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